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We make a Monte Carlo study ol the coupled two-scalar Xcj>icj>2 model in four dimensions at fi- 
nite temperature. We find no trace of Inverse Symmetry Breaking for values of the renormalized 
parameters for which perturbation theory predicts this phenomenon. 
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Is it always true that more heat means more disorder? 
The intuitive answer to this question would be yes and 
indeed this is what happens in the majority of physi- 
cal systems. Nevertheless, over 20 years ago, S. Wein- 
berg ||l|, quoting an unpublished remark by S. Coleman, 
observed that there may be exceptions to the general 
rule: in models with a sufficiently rich scalar sector, some 
of the scalars may acquire a negative Debye mass, with 
the result that the symmetric vacuum becomes neces- 
sarily unstable at high temperatures. This remarkable 
phenomenon is called Inverse-Symmetry-Breaking (ISB) 
or Symmetry-Non-Restoration (SNR) depending on the 
symmetry of the ground state at zero temperature. The 
possibility of ISB and SNR in realistic particle models 
and their cosmological consequences have been explored 
in a number of papers Weinberg's analysis of ISB 
and SNR was based on a simple one-loop approximation 
and some authors have questioned its reliability. Subse- 
quent studies, using different approximations with vari- 
ous amounts of non-perturbative content, produced con- 
tradictory results. While some concluded that ISB and 
SNR cannot occur others found that they do occur 
0, even though (in the majority of the cases) in a region 
of the parameter space significantly reduced with respect 
to the one-loop result . 

In this paper we present the first Monte Carlo study of 
ISB in 4 dimensions in a two-scalar model, with a global 
Z2 X Z2 symmetry. The job was mainly carried out on our 
RTNN computer, which holds 32 PentiumPro processors, 
for a total CPU time of approximately two months of 
the whole machine. The results were analogous to those 
found in 2+1 dimensions Q and seem to show that ISB 
is absent for certain values of the renormalized couplings 
in the region for which perturbation theory predicts it, 
in accordance with the theorem proven in Q . 

The model we simulated is described by the bare (eu- 
clidean) action: 
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If g^^'l < the potential is bounded from below for: 

>9.9(«)^ 5r^>0, 5f >0. (2) 

At high temperature, according to a simple one loop com- 
putation |l| , the thermal mass Mf (T) is proportional to 
{gf^ + g ) , where and are renormalized parame- 
ters. ISB (or SNR) occurs when the coupling among the 
two fields g^ is negative and such that, say. 
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Then, M|(r) becomes negative at high T and this means 
that if one starts from the disordered phase at T = 0, the 
field $2°'' should get ordered at sufficiently high temper- 
ature. This is the essence of ISB and what we have tried 
to investigate by means of a Monte Carlo simulation. 

As is customary on the lattice, we rewrite the action 
(|l|) in terms of dimensionless quantities as: 
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In principle, one would then proceed as follows. First, 
one would draw the phase diagram of the system (^ at 
T = 0. The physics described by renormalized pertur- 
bation theory would be recovered in the scaling region 
of this phase diagram, near the surface where both fields 
become critical. In order to check the prediction of ISB 
made by perturbation theory, one would then pick a point 
(the T = theory) in the disordered phase, well inside 
the scaling region. That point should be chosen such that 
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the renormalized couplings were as small as possible (for 
perturbation theory to be reliable) and verified the one- 
loop condition (^) for ISB as strongly as possible. For 
any choice of the bare parameters, the first condition is 
always satisfied if one goes deep enough in the scaling 
region, because the model is expected to be trivial. As 
for the second, it should be achieved by a suitable choice 
of the bare parameters. Having chosen the T = theory, 
one would then study the effects of temperature by sim- 
ulating the model on lattices having a finite extension Nt 
in the euclidean time direction. The physical tempera- 
ture is related to Nt as T — l/{aNt). For ISB to occur 
it should happen that the point for the T = theory 
selected above lied in the ordered region with (02) 7^ 
of all the phase diagrams with an Nt smaller than some 
Nt- On the contrary, the existence of ISB would be ruled 
out if one found that the size of the disordered region 
monotonically increases when the temperature is raised, 
namely when Nt is reduced, starting from Nt = oo, as it 
happens in normal systems, see for example 0|. 

In practice, it is impossible to realize this program, be- 
cause drawing five-dimensional phase diagrams with the 
high level of numerical precision that turned out to be 
necessary would have required an enormous time. Thus, 
we decided to fix once and for all the values of Ai, A2 
and A and then study the phase diagrams in the ki, K2 
plane. We took Ai = 0.3375, A2 = 0.01125. The condi- 
tion of stability for the bare potential, eq.(^), then gave 
a lower bound for A of —0.123, and for our simulations 
we selected the value A = —0.112 which is very close to 
the instability bound, thus giving us the best chances of 
observing ISB. 

The corresponding T = phase diagram, obtained on 
a single lattice, is shown in fig.|l|. We observe that 
the (ki,K2) plane is divided in four regions. The lower 
left region is the disordered phase, while the upper right 
one is the totally ordered phase. The wedges between 
them represent partially ordered phases: in the upper 
left wedge the field (j>2 is ordered, while (j>i is disordered, 
while in the lower right wedge it occurs the contrary. The 
critical line A is of second order, and we have obtained 
gaussian exponents. We have not studied accurately the 
thermodynamical limit of the B, D and E lines, but we 
saw a clear evidence of second order behavior on D and 
strong first order transitions across B and E. The critical 
lines cross at a point C, whose neighbourhood contains 
the scaling region because C evolves towards the gaussian 
fixed point when Ki — > 1/4 and Ai, A 0. Probing this 
region turned out to be practically impossible: while on 
small lattices it was not even possible to distinguish the 
various transition lines near C, on large ones the strong 
metastability due to the first order transitions prevented 
us from getting clean results. 

The next step was to check if in the scaling region of 
the disordered phase there were points where the renor- 
malized couplings satisfied the perturbative condition for 



ISB. We thus selected a point (ki = 0.255, K2 = 0.234) 
not far from the critical point C and measured there the 
renormalized couplings and the correlation lengths for 
both fields. 
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FIG. 1. Phase diagram at Ai = 0.3375, A2 = 0.01125, 
A = -0.112 for a 4* lattice. 

For these measurements we used Nf lattices, with 
Ns — 6, 8, 10, making between 50 and 250 million it- 
erations, with autocorrelation times of some tens of iter- 
ations. 

The renormalized couplings were measured using the 
following estimators: 
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where ^i{Ns, ki, K2) is the second moment correlation 
length 1^. Ui{Ns, Ki, K2) and Ui2{Ns, ki, K2) are the 
Binder cumulants: 
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where Mi is the average magnetization 
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1.570(2) 


2.784(2) 


32(4) 


2.6(7) 


-4.7(5) 



TABLE I. Correlation lengths and renormalized couplings 
at Ki = 0.255, ^2 = 0.234. 
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We checked the correctness of our measures by carry- 
ing out several tests in some weh-known models. In the 
first place we measured the renormalized coupling in the 
one-scalar 0'* theory, using the same lattices, and get- 
ting results fully compatible with those of ref . Q . In the 
second place we simulated the coupled system, but tak- 
ing (7^°^ = ^2"^ and 5^°^ = ffi^V^i which corresponds to 
the 0(2) theory, checking that these relations were also 
verified by the renormalized parameters. 

The results for our model are shown in table ^ The 
correlation lengths are stable with the lattice size, which 
indicates that we are using big enough lattices. On the 
other hand, they are sufficiently large to believe that we 
are in the scaling region. The renormalized couplings are 
rather stable as well. Though the large errors do not 
allow an accurate study of the thermodynamic limit, it 
is clear that eq.(^ is verified by g2 and g^, which is 
the essential thing. These renormalized quantities are 
all smaller (in modulus in the case of g^) than the bare 
ones {gf^ ~ 125, g^"^ ~ 4.9, g^°^ ~ -7.5), as expected 
from triviality, but while the values of (/|^ and g^ are 
rather small and well inside the perturbative region for 
ISB, gi is quite large. In order to get a significantly 
smaller value for gf^ we should have performed the mea- 
sures much deeper in the scaling region, but this would 
have been very difficult, because we would have needed 
much larger lattices and because, closer to the double 
critical point C, the nearby strong first order transition 
makes difficult to get clean results. 

The results of the numerical measurements were then 
extrapolated deeper in the scaling region using RG- 
equations, and taking the values of table | as initial 
data for the numerical integration. The relevant beta- 
functions were computed to two-loops order, keeping the 
corrections to scaling for finite values of the correlation 
lengths up to one-loop. In view of the large value of gf in 
the initial point, in the beta-function for g^, we included 
also the three loops contribution, of order ((jf^)'*. Accord- 
ing to 1^ , this approximation should be fully reliable for 
values of g^ even larger than ours. Upon integrating the 
RG equations along the line parallel to the ki axis pass- 
ing through the simulation point, we found that while gf" 
decreases rapidly, eventually entering in the perturbative 
region, the evolution of and is such that the per- 
turbative condition for ISB eq. (^ remains strongly satis- 
fied. We thus concluded that in the scaling region of our 
(ki,K2) plane there were points for which perturbation 
theory would have definitely predicted ISB for 02. 

Having analyzed the T — Q theory, we turned to the 
T > case. The issue was to study the direction of the 
shift of the critical line A as a function of the temper- 
ature. For that we needed an accurate determination 
of some critical points both on symmetric and on asym- 
metric lattices. Since they had to be found with a high 
numerical precision, in order to clearly distinguish the 



phase transitions corresponding to different values of Nt , 
we could not afford to explore the entire (ki, K2) plane, 
and thus we fixed once and for all the value of ki = 0.21, 
and searched on that vertical line for the accurate critical 
values of K2, Kl(A't), corresponding to the transition of 
(1)2. We will see that K^iNt) grows when Nt diminishes, 
i.e. that the disordered phase gets larger when increasing 
the temperature, indicating the non existence of ISB. 

Near the transition points, the correlation length for 
the field was about 0.7, not very large but still rea- 
sonable for us to believe that we were probing the scal- 
ing region. We simulated lattices with Nt = 2, 3, 4, 5, Ns 
and Ns = 6,8,10,12,14,16,20,24,28. More exactly, 
for Nt = 2,3,4,5,7V3, we simulated up to Ns"''''' = 
16, 20, 24, 28, 20 respectively, making sure in all cases that 
the results showed asymptotic behavior. 
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FIG. 2. Fits to obtain for Nt = 2, 5. 

For each of the points and lattices that were simu- 
lated we made between 5 and 15 million iterations of 
the Metropolis algorithm. The autocorrelation time T2 
ranged from 200 to 2,500 iterations. In some runs we 
used also the cluster algorithm, but found that it was 
not efficient for our values of the parameters. The errors 
in the estimation of the observables were calculated with 
the jackknife method. 

The estimators used and the method followed to take 
the thermodynamic limit were exactly the same as those 
used in ref. |^ and we refer to it the reader for the details. 
Having measured the Binder cumulants, we extrapolated 
them in a narrow K2 interval around the simulation point 
using the spectral density method. The values of k§ in 
the thermodynamic limit were obtained by looking at 
the intersections K2{Nsi, Ns2) among all possible pairs of 
curves Ui{Ns, K2) and using the following scaling law 

n*2iNs,bNs) - 4 = i^-^7V--i/- , (11) 

where lu is the exponent for the corrections to scaling. 
In fig.0 we show the fits to the above expression for the 
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Nt = 2, 5 lattices. For the calculation of the error, due 
account was taken of the fact that the pair crossings were 
not all independent of each other. 

For the asymmetric lattices, since the scaling parame- 
ter is Ns, while Nt is fixed, according to the hypothesis of 
dimensional reduction and universality, we used the ex- 
ponents of the Ising model in three dimensions, namely 
ly = 0.63 and uj ~ 0.8. As a check, we also computed ly 
directly from the data obtaining values fully compatible 
with the above one. 

For the symmetric lattices we should use the mean 
field exponents u! — 0, u ~ 1/2. In this case, apart from 
logarithmic corrections, all the crossings should occur for 
the same value of K2, and this is approximately what 
happens. In order to give an estimate of K2, we took 
the limit of the result for cj ^ 0, obtaining a very fast 
convergence. 

In this way we have obtained, for each value of Nt, an 
estimation of K2{Nt) with A^^ ^ 00. This constitutes the 
most important result of the simulations and is shown 
in figH: for a constant ki, starting from the symmetric 
lattices (the point with 1/Nt = 0), for the transition 
of the <j)2 field increases monotonically when Nt is de- 
creased; this means that the critical points for Nt finite 
shift deeper and deeper in the ordered region of the T = 
model, i.e. raising the temperature disorders the system, 
as it happens in normal cases (see for a comparison 
and contrary to what is required for ISB to occur. In 
order to make sure that the values of Nt that we simu- 
lated were large enough, we checked the scaling of 
function of Nt- From FSS analysis one would in fact 
expect a scaling of the form K2{Nt) — ^2(00) ~ Nt ■, 
with V = 1/2. Letting v variable, the best fit to the pre- 
vious expression is obtained for \/v = 1.8(2); this is the 
fit shown in the figure. This makes us confident about 
the fact that no changes in the behavior of {Nt ) should 
be seen for larger iVt's. 

In principle, our measurements are not conclusive, in 
that we cannot exclude a priori the possibility that the 
(/)2 transition lines of the finite Nt lattices move below 
that of the symmetric lattice for values of ki closer to C. 
Unfortunately we could not obtain accurate results in a 
reasonable computer time in that region, due to the ex- 
istence of strong metastabilities associated with the first 
order transitions. However, a toy-test on the 4^ and 4'^ x 2 
lattices near the point C confirmed the absence of ISB. 

In conclusion, we have found no trace of ISB in our 
two-scalar model, in a region of parameters for which 
perturbation theory predicts this phenomenon. 
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FIG. 3. K2 as a function of NjT 
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